Comment on "Spherical 2 + p spin-glass model: An analytically solvable model with 

a glass-to-glass transition" 
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Guided by old results on simple mode-coupling models displaying glass-glass transitions, we 
demonstrate, through a crude analysis of the solution with one step of replica symmetry break- 
ing (IRSB) derived by Crisanti and Leuzzi for the spherical s ~\-p mean-field spin glass [Phys. Rev. 
B 73, 014412 (2006)], that the phase behavior of these systems is not yet fully understood when s 
and p are well separated. First, there seems to be a possibility of glass-glass transition scenarios in 
these systems. Second, we find clear indications that the IRSB solution cannot be correct in the 
full glassy phase. Therefore, while the proposed analysis is clearly naive and probably inexact, it 
definitely calls for a reassessment of the physics of these systems, with the promise of potentially 
interesting new developments in the theory of disordered and complex systems. 
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In a recent paper^ , Crisanti and Leuzzi have proposed 
a detailed analysis of the so-caUed spherical s + p mean- 
field spin-glass model, defined by the Hamiltonian 
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where s and p are integers such that 2 < s < p, the spins 
CTi are A'^ real variables subject to the spherical constraint 



Sill = -^i ^'^'i random coupling constants J^"^ 
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and J, 



(p) 



are uncorrelated zero mean Gaussian vari- 



ables with variances slJ^/S^"^^ and p\j'^/2NP~^, re- 
spectively. At the inverse temperature [3, these variances 
enter into the definition of the control parameters of the 
model, /is — (/3Js)^s/2 and /ip — (/3Jp)^p/2. 

In fact, Ref. 1 mostly deals with the case s = 2, p > 4, 
which displays a rich phase diagram with four different 
phases and a variety of transitions between them^ . The 
case 2 < s < p, to which the present Comment is devoted, 
is only discussed shortly in an appendix and illustrated 
with the 3-1-4 model, whose static and dynamical phase 
diagrams both exhibit only two phases, paramagnetic 
with replica symmetry (RS) and glassy with one step of 
replica symmetry breaking (IRSB), separated by smooth 
transition lines. This result is claimed to be generic. 

In this Comment, we show that this might not be 
the case and that there could be a parameter domain 
apparently left unexplored by all previous workers and 
corresponding to 2 < s < p and p — s large enough, 
where new phenomena occur. To achieve this goal, we 
first recall seemingly little known results about simple 
schematic models of the mode-coupling theory (MCT) 
for the liquid-glass transition'^''*'^, which show glass-glass 
(G-G) transitions and higher-order singularities. Then, 
guided by the insight gained within the framework of the 
MCT, we report the results of a naive investigation of 
the IRSB solution derived by Crisanti and Leuzzi, which 
raise a number of issues about the structure of the phase 
diagram of the systems with large p — s. 

In the framework of the MCT-^'^'^, the schematic mod- 
els are minimal models reproducing the typical nonlin- 



earities and bifurcation scenarios of the dynamical equa- 
tions generically derived within this theory. A par- 
ticularly important class consists of the so-called Fpq 
models, 1 < p < q, defined by the memory kernels 
m{t) — Vp<f){ty + Vq<f){tY . They have been widely studied 
and, although most of the works have dealt with the F12 
and i^i3 models'^, scattered results are also available on 
more general Fpq models^^^. 

The bifurcation analysis of the solutions of the Fpq 
equation with 1 < p < q leads to a parametric represen- 
tation of the dynamical transition line, which uses the 
non-ergodicity parameter / = limt^+oo 0(^) as the vari- 
able. It reads 
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with {p—l)/p < f < {q—l)/q. Topological changes occur 
in this curve depending on p and q, or more precisely on 
their combination^ 
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If p and q are close enough, such that A{p,q) < 1, the 
curve {vp{f),Vq{f)} is smooth, but, if they are widely 
separated, such that A(p, 5) > 1, a loop appears in 
the glassy domain. Using then a dynamical stability 
criterion which states that, when several values of the 
non-ergodicity parameter seem possible at a given point 
(vp^Vq), only the largest one is physically meaningful, 
one shows that two of the three branches of the loop 
have to be discarded, while the remaining one forms a 
G-G transition line which terminates at an A3 higher- 
order singularity (the ordinary MCT bifurcation is of 
type A2)'^. In the marginal cases like the models with 
p = 2 and q = 9 ot p — 9 and q — 50, where A(p, q) = 1, 
the loop reduces to a point and yields an A4 singularity. 
These results are illustrated in Fig. 1 for different models 
with p ~ 2. Interestingly, exactly the same dynamical 
scenarios are obtained in MCT studies of colloidal sus- 
pensions with short-ranged attractions^'^". 

Let's now come back to the original spin-glass problem. 
One of the motivations for the study of the mean-field 
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FIG. 1: Phase diagrams of the Fpq schematic models with 
p = 2 and g = 7, 9, 11, 13, and 15 (from left to right). The 
solid lines are the true transition lines and the dashed lines 

arc the portions of the {'t'pi f). ''.'</(/)} parametric curves which 
arc excluded by the dyuaniical stability criterion. The higher- 
order singularities (A3 and A4) and the crossing points (C) 
delimiting the G-G transition lines are indicated by symbols. 



spherical spin-glass models with multispin interactions 
has been the finding that the time evolution of the spin 
correlation function in the Langevin dynamics of these 
systems^"'^'"'^^'"'^^'^'* is described at high temperature by 
schematic mode-coupling equations. Namely, the high- 
temperature dynamics of the s+p spin glass is described 
by the schematic model, with the replacement 

Vs-i — jJLs a-ud fp_i = /ip. This analogy turned out to 
be a fruitful source of new theoretical developments in 
the physics of structural glasses, allowing to complement 
the MCT, originally a purely dynamical theory derived 
under the assumption of equilibrium, with very inter- 
esting static^^ and aging^^ scenarios. So, since the Fpq 
schematic models with well separated indices show inter- 
esting and nontrivial bifurcation scenarios in their glassy 
domain, one might legitimately wonder whether similar 
features could be present in the glassy phase of the s+p 
spin glass as studied by Crisanti and Leuzzi^. 

We first consider the statics. Its method of solution 
within the IRSB scheme is described in detail in Ref. 1 
and we only recall the expression of the free energy per 
spin $ in the IRSB phase as a function of the overlap 
parameter q and of the usual IRSB parameter m. It reads 
—p^{q, m) = S00 + G{q, m), where Sqo is the entropy per 
spin at infinite temperature and G{q, m) is given by 

2G(a, m) = ^ + ^ + {m-l) ( ^q^ + !^qA + 
s p ^ ' \ s p J 

— ln(l — q + mq) H ln(l — q). 

m m 



In this expression, the coupling constants /ig and /7,p arc 
functions of q and m as well, obtained from the station- 
arity conditions dqG{q, m) = and dmG{q, m) = 0. 

In the same spirit as in the schematic MCT calcu- 
lation, one is led to consider the parametric surface 
{^J-s{q, m), ^p{q, m), G{q, to)}, which provides one with an 
implicit representation of the free energy as a function 
of the control parameters, and to look for topological 
changes when s and p are varied. As illustrated for the 
3-1-16 model in the insert of Fig. 2, where a part of this 
surface is plotted, such changes indeed occur when p— s is 
increased. A swallowtail, where the free energy as a func- 
tion of /U3 and jJie would be multivalued, is clearly visible 
above the surface. This looks very much like the loops 
met in the van der Waals theory of phase coexistence. 
Following this analogy, since the saddle-point calculation 
of the free energy requires that G should be minimized 
for given Hs and Hp, a crude way of dealing with this 
feature of the surface would be to discard it in order 
to define uniquely the free energy at a given state point. 
Then, this leaves one with a line of double points (a "bin- 
odal"), along which two distinct {q,m) pairs correspond 
to the same value of the free energy of the glassy phase. 
When this line is crossed, by following a path of constant 
/Us for instance, a discontinuity in {q,m) occurs, while 
the free energy remains continuous. At the level of the 
present discussion, this line appears as a candidate for a 
line of static G-G transitions, starting at a "triple" point, 
where the paramagnetic phase and two IRSB phases with 
TO = 1 and different values of q would have the same free 
energy, and ending at a "critical" point, where the dou- 
blet {q.,m) is such that dqiisdm^i-p — dmfJsdqfip = and 
dq{dqp.sdmP'p — dml^'s^q^Lp) = (tlic first cquafity char- 
acterizes "spinodals" corresponding to the edges of the 
swallowtail and meeting at the "critical" point). 

The resulting static phase diagram for the 3-1-16 model 
would then be the one outlined in Fig. 2. Its shape seems 
generic for well separated s and p, for instance, s = 3 and 
35 > 13, s = 4 and p > 23, or s = 10 and p > 105. The 
larger s, the larger p — s has to be for a candidate static 
G-G transition line to exist. 

We now turn to the IRSB solution for the dynamics 
of the s+p model. Following Crisanti and Leuzzi, the 
discussion is not based on a direct analysis of the off- 
equilibrium dynamics of the system, but on the study of 
its complexity E, a procedure which should yield identi- 
cal results. Thus it amounts here again to the considera- 
tion of a set of parametric equations, where ^s, ^p, and 
S are expressed as functions of the dynamical overlap 
and IRSB parameters. See Ref. 1 for the corresponding 
expressions. 

Like in the static calculation, the parametric surface 
{/Xs(g, to), /ip(g, to), i;(<7, m)} becomes singular for large 
p — s, with domains where 'E{p,s,fip) would be multival- 
ued. This is illustrated in Fig. 3 for the 3-1-16 model, 
where one can also see that the surface is more compli- 
cated than in the statics (at this point, it greatly helps to 
inspect the surface with an interactive computer graphic 
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FIG. 2: Putative IRSB static phase diagram of the 3 + 16 
model in the (/^Sj/iig) plane. T and C denote the "triple" 
and "critical' points and delimit the candidate G-G transition 
hne. Insert: Projection onto the {fi3,fj.ie) plane of the para- 
metric surface {jj,3{q,m), fiiQ{q,m),G{q,in)} in the region of 
the candidate G-G transition line. The outline of a swallow- 
tail is clearly visible. 



software in order to fully appreciate the nature of the 
singularities). Following again the analogy with the van 
der Waals theory, since the calculation of the dynamical 
solution requires that S should be maximized, a naive 
way to define a unique complexity at a given state point 
(/is, /ip) is to discard the parts of the parametric surface 
which do not correspond to the largest possible value of 
E. In the case of the 3-1-16 model, this reduction leads 
to the appearance of three distinct curve portions, suc- 
cessively delimited by points A, B, C, and D, in Fig. 3, 
which, when crossed, for instance along lines of constant 
/i3, are associated with discontinuities in {q,m). 

Line AB is a portion of the dynamical m — 1 line^^. 
Indeed, in the corresponding domain of /is and /xie, the 
sheet of the parametric surface originating in the plane 
/t3 = (the right part of the surface in Fig. 3) over- 
hangs the one originating in the plane /iie = (the left 
part). This means in particular that, when point A is 
approached along the border between the paramagnetic 
and glassy phases, the complexity obtained when coming 
from the large q side (from the right) is strictly larger 
than the one obtained when coming from the small q 
side (from the left). Thus, when the representative point 
of the system crosses the portion of the m = 1 line com- 
prised between points A and B, discontinuities in (g, m) 
(with m jumping to or from unity depending of the direc- 
tion of approach) and S occur, as the result of jumping 
from one sheet of the complexity surface to the other. At 
point B, the two sheets intersect for the first time. The 
m = 1 line continues below the reduced complexity sur- 
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FIG. 3: Top view of the parametric surface 
{fi3{q,m), fiie{q,m),'E{q,m)} for the dynamical IRSB 
phase of the 3-1-16 model. Except in the rightmost part of 
the figure where the surface is cut off in order to show an 
underlying swallowtail, the lower parts of the surface are 
hidden, so that only the point corresponding to the largest 
complexity is visible for a given (/i3,/ii6) pair. Points A, 
B, G, and D, marked by crosses on a white background, 
delimit the different branches of the putative dynamical G-G 
transition line (see text for details). 



face and becomes irrelevant (at least at the level of the 
present crude analysis). Point B marks the beginning of 
a line of double points, line BC, along which two distinct 
(q, m) pairs correspond to the same value of the complex- 
ity of the glassy phase. When this line is crossed, (q, m) 
is discontinuous, while S is continuous. Eventually, at 
point C, the line of double points terminates at a line 
of cusp singularities which extends towards the "critical" 
point D, providing a final line of {q,m) discontinuities, 
with a discontinuous S anew. 

As already suggested about the statics, these lines of 
(q, m) discontinuities appear as candidates for dynami- 
cal G-G transition lines. A simple condition on s and p 
for their existence is obtained by looking for stationary 
points in the dynamical m = 1 line, as they are neces- 
sary for the presence of a point like point A. Since the 
equations of this m = 1 line and of the transition line 
of the i^(s_i)(p_i) schematic model coincide, this condi- 
tion is just the same as in the MCT calculation, i.e., 
A(s — l,p — 1) > 1, with A defined above. Thus, candi- 
date dynamical G-G transition lines are expected in the 
systems with s = 3 and p > 11, s = 4 and p > 16, or 
s = 10 and p > 52, for instance. Generically, keeping s 
fixed and increasing p, the possibility of a G-G transition 
line first appears in the dynamics, then in the statics, so 
that there is a whole class of systems like the 3-1-11 model 
which only display a candidate dynamical G-G transition 
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line. Finally, there are still limiting cases like the 3 + 10 
and 10 + 51 models, where point A itself is a stationary 
point of the dynamical m — 1 line, to which the putative 
dynamical G-G transition line reduces. 

So far, we have only considered the topological features 
of the complexity surface in the dynamical IRSB solution 
and not the actual values of this function. It came as a 
surprise in the course of the present investigation that 
negative values of the complexity could be found in a 
large parameter domain. For a given s, this phenomenon 
is observed for p larger than the one needed for a singu- 
larity to appear in the complexity surface, so that it does 
not occur for all systems with a candidate dynamical G- 
G transition line. Considering again the 3 + 16 model, 
it is easily evidenced by computing E(g,m) along lines 
of constant m, including m = 1^''. It results that there 
is a whole region above the line joining points C and D 
in Fig. 3 and corresponding to rather small values of m, 
where the only available value of the IRSB complexity is 
negative (this is not the case on the m = 1 line which falls 
in the domain where T,{fis,lJ-p) is multivalued and where, 
for a given (/Zg, /Xp) pair, there is the possibility of finding 
a positive value of S on another sheet of the complex- 
ity surface). As pointed out by Crisanti and Leuzzi"'^'', 
this finding signals an unanticipated breakdown of the 
IRSB scheme and calls for a complete reexamination of 
the phase behavior of the systems with large p — s. 

Note that many of the above features of the dynamics 
of the s + p spin-glass models could have been antici- 
pated from the work of Caiazzo et al}^, who studied the 
off-equilibrium dynamics of a lattice gas generalization 
of the s+p spin-glass model. Indeed, for this richer and 
more complex system, these authors found the possibility 
of dynamical G-G transition lines and, along these lines, 
of dynamical singularities marking qualitative changes 
in the aging behavior of the system. They also mention 



that, deeper in the glassy phase, the "one-step replica 
symmetry breaking solution should not hold and, instead, 
a spin-glass-like aging dynamics should be found" . Ap- 
parently, and quite surprisingly, the possibility that the 
same could be true in the simpler s+p spin-glass model 
was not envisaged. 

In conclusion, guided by previous rc;sults on simple 
MCT models displaying G-G transitions and higher- 
order singularities^''', we have proposed a crude analysis 
of the IRSB solution derived by Crisanti and Lcuzzi in 
Ref. 1 which shows that the phase behavior of the spheri- 
cal s+p mean-field spin glasses with well separated s and 
p does not appear to be fully understood yet. First, there 
seems to be a possibility of G-G transition scenarios both 
in the statics and the dynamics of these systems. Second, 
there are clear indications that the IRSB solution can- 
not be correct in the full glassy phase^'^. Obviously, it 
results from the latter finding that the present analysis 
based on this IRSB solution, and also on probably naive 
analogies with the van der Waals theory of phase coex- 
istence, is necessarily at least partly inconsistent. How- 
ever, it is sufficient to raise a number of issues about the 
structure of the phase diagram of the systems with large 
p — s. It is beyond the scope of this comment and defi- 
nitely beyond the expertise of its author to answer these 
questions. So, we can only wish that the above observa- 
tions will be judged interesting by experts in the theory 
of spin glasses and will motivate future studies of these 
systems, with potentially interesting new developments, 
for instance, in the theory of the amorphous-amorphous 
transitions in structural glasses. 

As a non-expert in spin-glass theory, the author is in- 
debted to A. Crisanti and L. Leuzzi for useful advices on 
the presentation of this work and for insightful observa- 
tions on the physical implications of the reported results. 
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